Conditions for Phase Equilibrium in Supernovae, Proto-Neutron and Neutron Stars 
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We investigate the qualitative properties of phase transitions in a general way, if not the sin- 
gle particle numbers of the system but only some particular charges like e.g. baryon number are 
conserved. In addition to globally conserved charges we analyze the implications of locally con- 
served charge fractions, like e.g. local electric charge neutrality or locally fixed proton or lepton 
fractions. The conditions for phase equilibrium are derived and it is shown, that the properties of 
the phase transition do not depend on the locally conserved fractions. Finally, the general formalism 
is applied to the liquid-gas phase transition and the hadron-quark phase transition for typical astro- 
physical environments like in supernovae, proto-neutron or a neutron stars. We demonstrate that 
the Maxwell construction known from cold-deleptonized neutron star matter with two locally charge 
neutral phases requires modifications and further assumptions concerning the applicability for hot 
lepton-rich matter. All possible combinations of local and global conservation laws are analyzed, 
and the physical meaningful cases are identified. Several new kinds of mixed phases are presented, 
as e.g. a locally charge neutral mixed phase in proto-neutron stars which will disappear during the 
cooling and deleptonization of the proto-neutron star. 

PACS numbers: 26.50.+X 26.60.-c 05.70.Fh 64.60.-i 12.38. Mh 



I. INTRODUCTION 

Phase transitions (PTs) of the first kind play a crucial 
role in various astrophysical systems like cold neutron 
stars, proto-neutron stars, collisions of compact stars or 
supernovae. In all of these systems the equation of state 
(EoS) is one of the key information for the understanding 
of the various phenomena connected with them. For hy- 
drodynamic and hydrostatic calculations the EoS has to 
be known, and therefore the correct treatment of phase 
transitions is essential. 

In fact plenty of phase transitions can occur in the 
aforementioned systems: At densities below saturation 
density and temperatures lower than ~ 15 MeV the well- 
known liquid-gas phase transition of nuclear matter oc- 
curs, leading to a dense symmetric phase in coexistence 
with a more dilute neutron gas [l], HI, 0, 0, H| . At sev- 
eral times saturation density further phase transitions to 
exotic degrees of freedom like to quark matter in the chi- 
ral phase transition @, or to a kaon condensed phase 
[1, H, Qj| EI are possible. A phase transition to a pion 
condensed phase [12l fl3l . [l4| or to hyperon matter [l5[ 
might also be of first order. Phase transitions between 
different types of color superconducting quark matter 
were proposed in Refs. flU fl7l. [l8| . 

The inclusion of a phase transition to exotic degrees of 
freedoms can substantially alter the stability of the star. 
Usually, the phase transition leads to a softening of the 
EoS and therefore lowers the maximum mass which can 
be supported by the star. Thus, the theoretical predic- 
tions can be confronted with real observations of neu- 
tron star masses. Besides the mass, also other observ- 
ables can be linked to phase transitions in compact stars, 
e.g. gamma ray bursts [ljj, |20j, |2l|, |22j or sudden spin ups 
during the rotational evolution of young pulsars [23|, [24| . 
Furthermore the appearance and the structure of mixed 



phases can have important consequences for transport 
properties like the thermal conductivity or the neutrino 
cmissivitics and opacities [2^|. Also the shear modulus 
and the bulk viscosity will be altered, affecting the glitch 
phenomena or r- modes 0, [2?J . Consequently, the occur- 
rence of mixed phases can modify the thermal [28j and 
rotational evolution of compact stars. 

The description of phase transitions in cold, dclcp- 
tonized neutron stars is rather well understood and ex- 
tensively discussed in the literature. For bulk matter 
there are two, in principle different, possible treatments: 
the Maxwell construction based on local charge neutral- 
ity leading to a first order phase transition or the Gibbs 
conditions for phase equilibrium where only global charge 
neutrality is required. For the latter case, which was first 
discussed in [29(, the phase transition becomes second 
order according to Ehrenfest's classification and an ex- 
tended mixed phase appears inside the neutron star [3(| • 
If one wants to go beyond the bulk limit, finite size effects 
in form of surface and Coulomb e nerg ies need to be in- 
cluded, as was done in [?1 l3ll. l32l. l33l. 13-41] for a mixed phase 
of (hyperonic) hadronic matter and quark matter. Al- 
though Gibbs conditions are used in most of the publica- 
tions about the hadron-quark phase transition in the bulk 
approximation, as e.g. in Ref. [3^|, these detailed calcu- 
lations show that the EoS of the phase transition is in 
fact more similar to the Maxwellian case, if Coulomb and 
surface contributions arc included. Because of large sur- 
face tensions and small Debye screening lengths charge 
screening effects arc large, and the phases become al- 
most charge neutral. Then also global properties, like 
the mass-radius relation of cold neutron stars, resemble 
more the results of the Maxwell construction. Already 
in Q it was estimated, that for a > 70 MeV/fm 2 the 
Maxwellian case is recovered, recently validated by [3^ ]. 
In [36[ a first order phase transition to a kaon condensed 
phase was studied and similar results are found. For the 
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case of nuclear mixed phases (liquid-gas PT) the same 
charge screening effects are observed, but much less pro- 
nounced because of the low charge and baryon densities 
involved [37| . 

Recently, new works about the QCD phase transition 
in core-collapse supernovae or in compact star mergers 
appeared, as e.g. [H, HI] or [4fj]. For the proper descrip- 
tion of these more complex dynamical scenarios detailed 
numerical simulations are necessary. Within certain sim- 
plifications and/or model assumptions some interesting 
results were found. E.g. in Ref. [381 ] a successful explosion 
in spherical symmetry of a 15 M Q progenitor star was ob- 
tained due to an early appearance of a mixed phase of 
quarks and nucleons. 

For the early stages after the bounce in a supernova, 
one can represent the evolution of the proto-neutron star 
by some typical static configurations, see e.g. [U |42| . 
Until now, the conditions for phase equilibrium of matter 
in supernovae and proto-neutron stars were not discussed 
in detail in the literature. If only global charge neutral- 
ity is considered, the situation is clear and completely 
described in [2^]. However, in some cases the Gibbs con- 
ditions for phase equilibrium cannot be fulfilled at all, or 
only the simple Maxwell construction is wanted for the 
sake of simplicity, as e.g. in Ref. [43j], because then no 
mixed phase is present. But as will be shown in this arti- 
cle, the Maxwell construction has some subtleties and the 
usual procedure for cold deleptonized neutron stars can- 
not be used for proto-neutron stars. The requirement of 
conservation of lepton and/or proton number in addition 
to baryon number leads to significant differences in the 
equilibrium conditions, which was not taken into account 
in several previous publications, like e.g. in 0, S, [46[. 
The simple Maxwell construction (in the sense that the 
mixed phase vanishes) is only possible if in addition to lo- 
cal charge neutrality some other charges are fixed locally, 
resulting in particular new conditions for phase equilib- 
rium. Because of the additional conserved charges in- 
volved, there is a large variety of different descriptions 
of the phase transition, all of them representing different 
physical scenarios. We note that one of the new Maxwell 
constructions which we will present in this article, was al- 
ready considered in the work by Lugones and Benvenuto 

E3- 

The structure of this article is as follows: In section HT1 
we derive the equilibrium conditions for a general system 
of two phases in which not the single particle numbers, 
but only some global charges are conserved. In addition 
to globally conserved charges, we also consider the case 
in which some charge fractions are conserved locally. To 
achieve a general description which can be applied to all 
the combinations of global and local conservation laws, 
we start with a formulation, in which all the charge frac- 
tions are fixed locally. Then we derive the new equi- 
librium conditions if some of these local constraints are 
lifted. In section IIIII we present the properties of the 
phase transition depending on whether one has a singlc- 
or a multi-component body (only one or multiple glob- 



ally conserved charges). While most of this was already 
discussed in [2^] we want to reformulate it by using con- 
sistently our formalism and notations. Quite remarkably, 
we will show that locally conserved fractions do not in- 
fluence the qualitative behavior of the mixed phase. In 
section |IV] the general results will be applied to single and 
mixed phases of hadronic matter and/or quark matter. 
We differentiate between isothermal and adiabatic phase 
transitions. All relevant possibilities of local and global 
conservation of the different conserved charges will be 
analyzed and the corresponding equilibrium conditions 
will be presented. We will use the results of the previous 
sections and will show which assumptions are necessary 
to achieve the simple Maxwell construction even in the 
case of multiple conserved charges. Some of the combi- 
nations of local and global conservation laws have not 
been discussed in the literature so far, leading to new 
conditions for chemical equilibrium and new interesting 
properties of the mixed phase. In addition we discuss 
the role of neutrinos: We will show explicitly that also 
in a mixed phase the neutrinos can be handled as an 
independent contribution to the EoS, as long as the pro- 
ton fraction is conserved. At the end we will discuss the 
role of the additional quantum number of strangeness in 
strange matter, i.e. three-flavor quark matter or hyper- 
onic matter. In section [V] we will summarize the main 
results and draw the conclusions. 

II. EQUILIBRIUM CONDITIONS 

Consider a thermodynamic system with volume V and 
temperature T composed of two different phases. The 
numbers of the N 1 (N 11 ) different particle species of 
phase / (//) are denoted by N? {Nj 1 ). To distinguish 
the two phases we introduce the index k = I, II. The 
thermodynamic potential of the system is the Helmholtz 
free energy F(T, V, Nf). 

In many cases, there exist some conserved charges, but 
the single particle numbers are not conserved. Then it 
will be more convenient to use these conserved charges as 
the independent degrees of freedom for the description of 
the state of the system, instead of specifying all the single 
particle numbers. In this section the equilibrium condi- 
tions in terms of the chemical potentials of the particles 
shall be derived. 

Let us assume there are C conserved charges C a , a = 
0, C — 1 for which a conservation law of the following 
form exists: 

C^N^Nj 1 ) = Ci{Nl) + C^(7Vf ) = const. , 

C£(Atf) = E Q »^ (!) 

i 

with a\ denoting the amount of charge C„ carried by 
particle i of phase k. The total conserved charge consists 
of the charge in phase I and in phase II. 

In addition to the C global conservation laws, we con- 
sider additional local constraints, which depend only on 
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the particles in one of the two phases. To be more 
specific, we require that some (denoted by the index r, 
t ^ 0) of the local charge fractions of the globally con- 
served charges C a have to be equal in the two phases: 



Yj = K 



it 



Y k (N k ) = § = ^ £. a k T N k = const. 



(2) 
(3) 



we find: 



Y T = C T /Co = const. 



Y k = Y T = const. 



(5) 



Eq. §5$) is the reason why we can evaluate the local con- 
straint of cq. @ independently for the two phases. 

To achieve a general description applicable for all kind 
of local constraints we will first use the two local parts Cq 
of Cq and all the local charge fractions Y k , t = 1, C — 
1, as the state variables, i.e. we treat them as indepen- 
dent degrees of freedom. In the following we want to 
derive the equilibrium conditions for the internal depen- 
dent degrees of freedom N k if only these state variables 
are specified F(T, V, N?) = F(T, V, N k (T, V, C l , Y}) and 
are kept constant. Thus: 



Yi 



const, 
const. 



(6) 
(7) 



Later we connect these state variables with the real local 
constraints, in which only some of the local fractions have 
to be equal and Cq is conserved only globally. 

From the first and second law of thermodynamics we 
get for the free energy F expressed by the particle num- 
bers Nf 



= dF 



E 



dF 



dN k 



dm 



(8) 



if the volume and the temperature are kept constant, k 
denotes the phase different to k. As the temperature 
is one of the state variables of the two phases it is set 
equal by construction, so that thermal equilibrium be- 
tween the two phases is assured. Only the total volume 
V = V 1 + V 11 is kept constant, but the two subvolumcs 
can vary, leading to pressure equality as the condition for 
mechanical equilibrium: 



P 



P 



ii 



With 



dN k 



F(N l j ) 



(9) 



(10) 



cq. (|5J) becomes 



i . k 



(ii) 



The constraints of eqs. © & ([7]) can be implemented by 
the means of Lagrange multipliers Aq , Ay- , by adding 



where Cq is the local part of the conserved charge Co and 

C k the local part of C Tl which is also one of the globally an( j 

conserved charges. Cq shall be a positive, non-vanishing 

quantity so that it is suitable to characterize the size of 

the phases in the mixture. If we introduce the global 

fraction Y T , 



^dC^=X k J2iaidN^ = 



\^dY k = A^g^dC* - Y k dC k ) = 



oLdNf 



(12) 



0, (13) 



(4) to dF, where we used eq. ([5]). This leads to: 



,,k \k„k | \k 1 (k v „k \ 

°0 



(14) 



dF can also be expressed as a function of (Cq, Y k ): 
dF 



Q = dF = £ 

k 

E 



dC k Q 
dF 



dC k 



dY k 



r<l vk -vk 



dYt 



(15) 



We introduce the chemical potential of the conserved 
charges Cq, 



d 



dC k 



fk-yl 



F(C l X) 



and of the conserved fractions Y k : 
d 



BY* 



(~<l vfc V fc 



F(C X) 



(16) 



(17) 



Then it is easy to realize, that the Lagrange multipliers 
are equal to the chemical potentials of the corresponding 
charges: Aq = yUg, A y = [i\ , so that: 



Y T a k ) 



(18) 



It is interesting to see, that the chemical potentials 
of the particles depend now directly on the value of the 
locally fixed charge fractions and the unknown value of 
Cq. 1/Cq appears because a change in Y k implies a 
change in the particle numbers proportional to Cq . With 

dF 



/4 



C, 



dC k 

(~<k k 



r<l vfc V fc 



(19) 



one can see that (Jy is proportional to Cq and the chem- 
ical potential of the charge C k . In fact the local chemical 
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potentials of the particles can only depend on other local 
intensive variables, so Cq has to drop out of eq. (fT5|) . 

By using this result we get the following expression for 
the equilibrium conditions for the chemical potentials of 
the particles: 



M; 



Em? 



(20) 



The first two of the three terms simply state that the 
chemical potential of particle (i, k) is given by the sum 
over the amount of conserved charges that the particle 
carries multiplied by the corresponding chemical poten- 
tials. The term proportional to Y T appears only for par- 
ticles which contribute to Cq. It is due to the change 
implied by dCp in the charge C k if Y k is kept con- 
stant. It would not appear if instead of the fractions 
the chargeswere used for the description of the state of 
the system. This shows the importance in the definition 
of the chemical potentials of which other quantities are 
kept constant. 

It is important to realize that the local chemical poten- 
tials Moi Mr wn l m general be different in the two phases, 
leading to different chemical potentials of all particles. 
Only particles of the same phase which carry the same 
quantum numbers will have equal chemical potentials. 

All together there are 2C+N 1 +N 11 + 1 unknown vari- 
ables in eq. ([20]) : the chemical potentials fj,Q , ii k , /j, k and 
one of the two subvolumes V . They can be determined 
from the N 1 + N 11 chemical equilibrium conditions (f2"0")) . 
pressure equilibrium ([9]) and the 2C conservation laws 
© and for the fixed state (Cfi,Yf). If all the re- 
lations N k = N k {T,V k , fij) are known, the system is 
determined completely. 

Eq. ([20]) can also be applied for a single phase k, by 
setting V k = for the other phase, equivalent to N k = 
0. Then pressure equilibrium is not required any more, 
and the whole system of equations can be solved and 
all thermodynamic quantities can be determined, too. 
If the number of conserved charges C is equal to the 
number of particles N k in this phase, then the conserved 
charges directly fix all the N k particle numbers N k . If 
C < N , N k — C equilibrium conditions between the 
chemical potentials of the particles will exist. 

Next we want to understand the consequences if in- 
deed only Cp but not the Cp are conserved and if be- 
sides the global conservation of the other charges C T , 
r = 1, C — 1 only for the L fractions Y k local con- 
straints of the form exist. I.e. we want to connect 
the chosen state variables to some particular set of local 
conservation laws. We will denote the only globally con- 
served charge fractions by Yy, j = 1, G — 1, with the 
number of globally conserved charges G. For the index 
A of the locally conserved charge fractions Y k we then 
have A = G, C - 1 so that C = G + L. 

The global constraint for Cp can be written as: 



ApdCo = A > ; a k dN k = 



i , k 



(21) 



All the global charge fractions Y 1 
served: 



C 7 /Cp are also con- 



Ay, ^-(dC 7 - y 7 dC ) = 



t Co 

a 



k dN k 



Y^i^dNf =0.(22) 



We already implemented the new Lagrange multipliers 
Ao, Ay 7 ■ 

By comparing eqs. (H2J) & (JTSJ) with (JUJ) fc J2"2"|) . one 
finds that the local conservation laws lead to the same 
constraints (|2"Tj) & (|2"2"|) which are added to dF if we set: 



Aq 



A - 


C' 



A 

CP 



An 



Co 



which is equivalent to: 



Mo — Mo J — : Mo 



c' n ' 



«*■ M 7 



M 7 



(23) 
(24) 

(25) 
(26) 



where we used eq. (TK)|) in the last line. Eqs. (|2"5|) & 
(|26p are the new additional equilibrium conditions for the 
globally conserved charges in terms of the local chemical 
potentials. To express the equality of these chemical po- 
tentials we introduced the variables /ip, /i 7 , which are the 
global chemical potentials of the corresponding charges. 
If a local conservation law is lifted, an additional con- 
dition for chemical equilibrium between the two phases 
appears, and the whole set of equilibrium equations can 
still be solved. Eqs. (f25|) & ([26]) are the expected result 
that local chemical potentials become equal in the two 
phases when the corresponding local constraint is lifted. 
Then the two local fractions can adjust to minimize the 
free energy. After equilibrium is reached, the two phases 
can be separated from each other without any changes 
arising. 

If one of the globally conserved charges, denoted by C$ 
in the following, is actually not conserved any more this 
means that 



dF 
do's 



= 



to minimize the free energy, leading to 

M5 = /4 7 = o . 



(27) 



(28) 



A non-conserved fraction gives two local constraints for 
the chemical potentials. The two eqs. (|28p replace the 
equilibrium condition (|2l)|) . With this additional infor- 
mation the whole system can be determined, even with- 
out knowing the value of C$ any more. This is in ac- 
cordance with our previous conclusion that all thermo- 
dynamic quantities can be determined, independently of 
the number of conserved charges C and the number of 
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particles N 1 and N 11 . We note that with the new in- 
formation of eq. I|28p the chemical potentials of the re- 
maining conserved charges can possibly be written in a 
different simplified form. 

This procedure for non-conserved charges can also be 
applied for a single phase k, in which only one chemical 
potential /j, k exists, leading to: 

$ = . (29) 

All other conclusions are also analog to the mixed phase. 

From cqs. (f20| it follows immediately, that two parti- 
cles of the same phase have equal chemical potentials if 
they carry the same quantum numbers: 

a4 - m« if i = &L • (so) 

If particles m and n of two different phases carry the 
same quantum numbers, = a„„, e.g. if some of the 
particles in the two different phases are identical, this is 
no longer true in general. The local chemical potentials 
of local fractions will be different in the two phases. Con- 
sequently, if a particle carries global and local charges, its 
chemical potential will also be different in the two phases. 
Only if they do not contribute to the locally conserved 
charges from eqs. flU), (E9 & (23) it follows that the 
chemical potentials of such particles are equal: 

\L = ■ (31) 

This means that, since such particles can be exchanged 
freely between the two phases, in equilibrium always the 
same amount of energy is needed when the number of 
particles N k is varied in one of the two phases. If no 
local constraints are applied, the chemical potentials of 
all identical particles become equal, which are the well- 
known Gibbs conditions for phase equilibrium. 

For fixed temperature T, pressure p and particle num- 
bers N k , the correct thermodynamic potential is the 
Gibbs potential G, which is also called the Gibbs free 
enthalpy. With eq. (HI) & (2H]) we get the following re- 
lation inside the mixed phase: 

G(T,p,N?(C ,Y T )) = J2^ N i 

i . k 

= Mo Co = ^C / + ^ / C7 /7 ,(32) 

independently of which local constraints are applied. For 
a single phase k with Cq = C and V k — V, eq. PD|) 
leads to: 

G(T,p,N?(C Q ,Y T )) = = ^Co ■ (33) 

i 

In this case in principle the index k can also be suppressed 
because only one single phase exists. These two relations 
can also be used in other thermodynamic potentials or 
in the fundamental relation of thermodynamics. 



III. PROPERTIES OF THE PHASE 
TRANSITION 

Depending on the number of globally conserved 
charges the properties of a phase transition are quali- 
tatively different. The differences between a single and 
a multi component body were extensively discussed in 
f29| | for neutron stars. In the following we will show that 
locally fixed charge fractions do not influence the quali- 
tative behavior of the phase transition. 

For numerical simulations of proto-ncutron stars and 
supcrnovae, usually an equation of state in tabular form, 
as e.g. gl El S3], in terms of (T,co,Y T ), r = 1,...,C-1, 
with the density cq = Co/V is applied (in these EoSs the 
state variables are the temperature, the baryon density 
and the proton fraction Y p ). If the volume V (which 
is completely arbitrary in the thermodynamic limit) is 
also taken as one of the state variables, all the conserved 
charges C a , a — 0, ...,C — 1, are known, too, and the 
Helmholtz free energy becomes the appropriate poten- 
tial: F = F(T, V, C„). In the following we want to dis- 
cuss the properties of the mixed phase for such an EoS 
in terms of (T, V, C a ) if only V is varied and the other 
state variables C a and T are kept constant, i.e. we are 
investigating an isothermal change of cq at constant Y T . 
As before we assume that only L < C — 1 of the frac- 
tions are fixed locally in the form: Yl = Y^ 1 = Y\. For 
the globally conserved fractions we will continue to use 
the index 7 instead. The number of globally conserved 
charges/fractions is then G = C — L. 

By analyzing further the condition for pressure equi- 
librium one can make statements about the qualitative 
properties of the phase transition and the possible ap- 
pearance of mixed phases. The pressure in each phase 
can only depend on the local chemical potentials of the 
particles in this phase: 

p I (T,fxl)=p II (T,^ j I ). (34) 

We already argued before, that also for a single phase 
the knowledge of the conserved charges, the temperature 
and the volume V is sufficient to determine all thermo- 
dynamic quantities. Obviously, the chemical potentials 
of the particles cannot depend on the size of the phase, 
so that they have to be determinable by the local density 
Cq and the local fractions Yf alone: 

p k (T,tf)= P k (T,c%,Y T k ). (35) 

Next, we can use [1q instead of c k , and replace the un- 
known local fractions of the fractions Yy which are con- 
served only globally by their local chemical potentials: 

p k = P k (T,[i k ,tx k ,Y£). (36) 

According to cqs. ([25)) & (|26|) the remaining chemical 
potentials have to be equal in the two phases, and the 
locally fixed fractions, too. Thus: 

;/(T,^ ,/i 7 , Y~ A ) =p II (T,no,fj, 7 ,Y x ) . (37) 
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This formulation, in which the chemical potentials of the 
globally conserved charges are treated as known state 
variables, is most convenient to discuss the properties of 
the phase transition. Of course, the chemical potentials 
in this equation are actually fixed by the values of the 
density Co, the fractions Y T and the chosen local con- 
straints. 

First we will analyze the case in which no other globally 
conserved charges besides Co exist, G = 1. Eq. (pTT|) then 
leads to a relation 

f i =ftf ex (T,Y x ), (38) 

which also fixes the coexistence pressure: 

p = p coex (T,Y x ) . (39) 

This means that there is only one value of the pressure 
pcoex anc j £ ne cnem i ca j potential /ig oe:E , where the two 
phases can coexist. All other local intensive variables are 
also fixed by the local constraints and the equilibrium 
conditions and remain constant in the mixed phase, too. 
A simple way to determine the phase transition pressure 
is to see where the p k (T, //o, YA)-curves of the two phases 
intersect. The transition from phase I to phase II occurs, 
when the pressure is equal in the two phases. 

The mixed phase will extend over a certain range in the 
density Co, with the onset given by Cq(T, ^Q° ex , Y T ) and 
the end by Cq / (T, fJ,Q° ex , Y\). Inside the mixed phase, the 
intensive variables are independent of Co and Co is only 
used to specify the volume fraction < X = V 11 /V 1 < 1 
of the two phases: 

c = (1 - X)cl{T, Y X ) + X<g(T, nT\ Yx) ■ (40) 

All extensive variables change linearly with the volume 
fraction in the same way, e.g.: 

e = (1 - X)4(T, ^° eX ,Yx) + Xe^iT, ^,Y X ) . (41) 

Therefore the calculation of the mixed phase becomes 
trivial. After the coexistence pressure is found it is given 
by a linear interpolation in the volume fraction between 
the onset and endpoint of the mixed phase. This case 
corresponds to the well known Maxwell construction. 

For G > 2 globally conserved charges, cq. (|37[) is not 
sufficient to determine all chemical potentials |«o, £t 7 - The 
equilibrium condition only allows to fix one of the chemi- 
cal potentials, e.g. /Uq = fJ-o(T, Y T , /x 7 ) and G — 1 chemical 
potentials remain unknown. But besides the equilibrium 
conditions, also the total volume and the globally con- 
served charges have to have the correct value: 

V = V 1 + v n 
C 1 = Cl / (V I ,T,n ,^,Y x ) + 

^ 7 (F /7 ,T> ,M7^a) ■ (42) 

These G + 1 cqs. involve only two further unknowns 
V 1 and V 11 so that the whole system of eqs. (|3T|) and 



(|42p can be solved for given volume V, and all thermo- 
dynamic variables can be determined. Consequently in 
this case all quantities (including the pressure) will de- 
pend on the values of the densities Co and c 7 = C 7 /V. 
A change in the density cq will also imply a change in 
the pressure. Thus for G > 2 there will be an extended 
range in pressure in which the two phases can coexist. 
The simple Maxwell construction cannot be applied, as 
the system does not behave linearly any more. Instead it 
is necessary to calculate the mixed phase at every point 
(T, V, Co) explicitly. 

In both cases (G = 1 or G > 2) an extended mixed 
phase between the two phases forms for the chosen state 
variables. At the onset of the mixed phase the volume 
of the newly appearing phase vanishes. Similarly, at the 
end of the mixed phase only the second phase remains 
(we assume that the two EoSs of the single phases do 
not show any discontinuities). Thus the mixed phase be- 
comes identical to the neighboring single phases when ap- 
proaching the onset or end of the mixed phase. Inside the 
mixed phase the volume fraction changes continuously 
from to 1, and consequently all global thermodynamic 
variables up to first derivatives of the thermodynamic po- 
tential will change continuously across the whole phase 
transition. The second derivatives will in general be dis- 
continuous, as they involve the derivative of the volume 
fraction. 

The thermodynamic potential is the Hclmholtz free en- 
ergy, which also changes continuously and which has the 
following form inside the mixed phase: 

F = -pV + ^NM 

i,k 

= -pV + MoO), (43) 

where we used eq. (|32|) . For G = 1 in which the pressure 
and the chemical potentials arc constant, the free energy 
changes linearly with the volume V. 

Next we want to use the different set of state variables 
(T, p, C a ), in which the volume is replaced by the pres- 
sure and we want to analyze the properties of the phase 
transition in this new formulation. Now the Gibbs free 
enthalpy G = G(T,p 7 C a ), already specified by eq. (|3"2"j) . 
is the appropriate thermodynamic potential. These state 
variables are especially important because they can di- 
rectly be used for the description of isothermal neutron 
stars. Under the influence of gravity in a hydrostatic con- 
figuration of a compact star, the pressure has to change 
continuously and has to be strictly monotonic. In the 
following we only consider a change of the pressure p and 
keep all the other state variables constant. 

If G = 1 the mixed phase collapses to one single 
point at the coexistence pressure p coex introduced above. 
There is only a point of coexistence, but no extended 
mixed phase. No mixed phase has to be calculated, only 
the transition point has to be determined. As /io is con- 
stant across the mixed phase and continuous at the end- 
points, one gets that the potentials of the two phases are 
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equal at the transition point, G = G 1 = G 11 . The equal- 
ity of //q = /Xq 7 leads to the equality of the Gibbs free 
enthalpy. 

This equality can also be seen as the reason why the 
volume fraction of the two phases remains arbitrary at 
the coexistence point and cannot be determined from the 
equilibrium conditions. Thus all extensive quantities but 
those of the externally fixed state variables remain un- 
specified at the coexistence point. 

Because of mechanical, thermal and (at least par- 
tial) chemical equilibrium the thermodynamic potential 
changes continuously across the transition, even though 
no extended mixed phase exists. For smaller or larger 
pressures the mixed phase disappears, and only the phase 
with the lower Gibbs free enthalpy will be present. De- 
spite this, the phase transition is not continuous, as 
e.g. the volume behaves discontinuously due to the dis- 
appearance of the mixed phase: 



lim 



dG 1 



dp 



V 1 ± V 



n 



lim 

p> >pCO 



dG 



n 



op 



T.,C a 



(44) 

Therefore the charge densities, defined by C a jV will 
change discontinuously, too. Also the entropy jumps in 
an analogous way at the phase transition, if T ^ 0. The 
internal energy, given by E = G — pV + TS will also 
behave discontinuously in general. These discontinuities 
appear in the first derivatives of the thermodynamic po- 
tential. Thus the phase transition is of first order ac- 
cording to the Ehrenfest classification. The discussed 
scenario is the familiar case of the Maxwell phase tran- 
sition, e.g. known from the liquid-gas phase transition of 
water. 

For G > 2 there will be an extended range in pressure 
in which the two phases can coexist and an extended 
mixed phase forms. As noted before, the simple Maxwell 
construction cannot be applied. The mixed phase does 
not behave linearly any more and thus it has to be calcu- 
lated explicitly for every single pressure. Now the equilib- 
rium conditions and the knowledge of the state variables 
become sufficient to specify the volume fractions and all 
other thermodynamic variables of the two phases. This 
case is usually called the Gibbs construction in the con- 
text of cold dclcptonizcd neutron stars with global charge 
neutrality. The presence of a mixed phases with X = at 
the onset and X = 1 at the endpoint assures that all ther- 
modynamic variables (up to first derivatives) change con- 
tinuously across the phase transition, as argued above. 
The phase transition is of second order according to the 
Ehrenfest classification. 

As was shown, locally conserved charge fractions do 
not influence the qualitative behavior of the phase tran- 
sition. It is only the number of globally conserved charges 
which determines the order of the phase transition. In- 
dependently of any locally conserved charge fractions, 
for G = 1 the system behaves like a simple body. By 
replacing globally conserved charges by adequate local 
conservation laws, this allows to reduce the number G 



of globally conserved charges. In the next section wc 
will use this procedure for isothermal PTs to arrive at a 
Maxwell construction, even if in principle multiple con- 
served charges exist. 

One can expect that the extension of the mixed phase 
decreases with the number of local constraints applied. 
When Co remains as the only globally conserved charge, 
the mixed phase will disappear completely inside a com- 
pact star. Furthermore the discontinuity of the second 
derivatives will increase with the number of local con- 
straints. For G = 1 even the first derivatives become 
discontinuous then. 

The situation becomes different, if we consider a phase 
transition in an adiabatic instead of an isothermal pro- 
cess. To keep the entropy constant, necessarily the tem- 
perature has to change across the transition with equal 
temperatures in the two phases at each point of the mixed 
phase. The change in temperature will lead to a change 
in all other intensive variables, too. Thus even for G = 1 
the Maxwell construction cannot be applied. If mechan- 
ical and thermal equilibrium between the two phases is 
required, the mixed phase docs not vanish in a hydro- 
static situation under the influence of gravity, even if only 
one globally conserved charge exists. This is an impor- 
tant result, especially for proto-neutron stars which are 
very often described by a constant ratio of entropy per 
baryon as a first approximation. Some other local con- 
straints are needed (e.g. locally fixed entropy per baryon) 
which will influence the conditions for phase equilibrium 
in a non-trivial way. Therefore wc will discuss adiabatic 
phase transition at the end of this article separately 



IV. APPLICATION FOR MATTER IN 
SUPERNOVAE, PROTO-NEUTRON STARS AND 
NEUTRON STARS 

As an example, the general relations which were found 
shall now be applied to the liquid-gas phase transition 
of nuclear matter and the hadron-quark phase transi- 
tion under typical astrophysical conditions. The hadronic 
phases shall consist of N y neutrinos, N e electrons, N p 
protons and N n neutrons (net numbers, including an- 
tiparticlcs). The two-flavor quark phase shall be com- 
posed out of iVe electrons, N v neutrinos, N u up and Nd 
down quarks. Furthermore, at the end of the section 
we will discuss strange quark matter, too, in which N s 
strange quarks are also part of the system. 

The baryon number Nb = N n + N p for hadrons, and 
Nb = 1/3(N U + Nd) for quarks, and the total electric 
charge number N c = N p - N e , and N c = 1/3(2N U - 
Nd)—N e , always have to be conserved in a closed system. 
Because of charge neutrality Nc = 0, but the concrete 
values of the conserved charges are actually irrelevant 
for the equilibrium conditions of the chemical potentials. 
Furthermore, there are two additional conserved charges 
possible, the lepton number Nb = N e + N v and the pro- 
ton number N p . The latter in combination with baryon 
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number conservation is equivalent to the conservation of 
isospin or of baryonic electric charge. Obviously, this is 
equivalent to flavor conservation for up and down quark 
matter. Thus, constant N p = 1/3(2N U — Nd) can be used 
to express isospin conservation for both kind of phases. 
To achieve a consistent description of quark and hadronic 
matter we will use the term proton number or proton 
fraction conservation as a synonym for the conservation 
of isospin in the two-flavor quark phase. 

Usually instead of fixing ( Nb , Nc , Nl ,N p ), an inten- 
sive formulation in terms of the proton and lepton frac- 
tions Y p = n p /riB and Yl = (n e + n v )/nB, the baryon 
number density n b and charge density nc = n p — n e = 
together with a fixed temperature T is used, like e.g. in 
Refs. [H, 0, We will also apply this formulation in 
the following. Instead of the electric charge per baryon 
Yc = used before we can also take nc — to describe 
the state of the matter. In the thermodynamic limit, the 
size of the system becomes irrelevant, so that we can as- 
sume that the volume V is also known. Obviously then 
it is completely equivalent to fix (T, ng, nc, Y p , Yl, V) 
instead of (T, N B , N c , N p , N L , V). 

Usually for fixed Y p the neutrinos are not included in 
the construction of the EoS but are treated separately. 
The neutrino dynamics play a crucial role in supcrnovae 
and proto-neutron stars. To describe the evolution of 
such systems it is necessary to handle the neutrinos with 
a detailed dynamical transport scheme in which their 
emission, scattering and absorption is calculated. We will 
explain why it is possible to construct the non-neutrino 
EoS without taking the neutrinos into account, by show- 
ing explicitly that the non-neutrino EoS is independent 
of the neutrino contribution in single phases and in mixed 
phases, too. 



A. Single Homogeneous Phases 

In table U the chemical potentials of the conserved 
charges Nb, Nc, N p and Nl are expressed in terms of 
the chemical potentials of the particles of one phase if 
the baryon number and all fractions are kept constant 
locally, leading to the unusual form of n k NB . As long as 
only one single homogeneous phase exists, local conser- 
vation laws are identical to global ones. The index k can 
be suppressed if we want to use table U for one single 
homogeneous phase. 

The special form of the chemical potentials of the con- 
served charges/fractions can be understood in a simple 
way. fic^ gives the change of the free energy with the 
change of the charge C 1 for constant proton and lepton 
fraction and electric charge neutrality. The combination 
of chemical potentials of the particles which is found for 
uc y corresponds to the change of the particle numbers in- 
duced by the change of C 7 under the chosen constraints. 
We note once more, that the form of a chemical poten- 
tial depends on which other quantities arc kept constant. 
E.g. the baryon chemical potential would be equal 



to /i„ (for nuclear matter) if instead of the fractions the 
charges themselves were used as the state variables. How- 
ever, the final equilibrium conditions are not (and cannot 
be) affected by the choice of the state variables. Thus we 
can use the description presented here, which is most 
convenient for our purpose as it can be applied for single 
phases as well as for all possible combinations of locally 
conserved fractions inside mixed phases. 

If all the four charges are conserved this corresponds 
to the situation of matter with completely trapped neu- 
trinos, but too short dynamical timescales to change the 
proton number by weak reactions. Next we will discuss 
the different possibilities in which some of the charges are 
actually not conserved any more. As was shown before, 
for every charge becoming not conserved an additional 
equilibrium condition appears. With the new informa- 
tion of eq. then the chemical potentials of the re- 
maining conserved charges can possibly be written in a 
different simplified form. We note that baryon number 
and electric charge always have to be conserved. 

For non-conserved lepton number from eq. (|28p and 
table Q] \i v = follows. Neutrinos are completely un- 
trapped/frce streaming. If the proton number is still con- 
served weak reactions involving nucleons are assumed to 
be completely suppressed. We will discuss the meaning 
of the result a v = in more detail later. 

If neutrinos are completely trapped and there is enough 
time for the weak reactions to become efficient, the lepton 
number is conserved but the proton number not. Then 
the well known weak equilibrium conditions 

- - flu + A*e = (45) 
for nuclear matter and 

- - + Me = (46) 

for quark matter are found. 

If then at a later stage in the evolution the neutrinos 
become completely untrapped, the lepton number is not 
conserved any more. Only baryon number and electric 
charge remain as conserved charges. Then two equilib- 
rium conditions are necessary to derive all particle num- 
bers Ni. Without lepton number conservation \i v = 
and the neutrinos drop out in the /3-equilibrium condi- 
tions 

u e + fi p - fi n = (47) 

for nucleons and 

u u - Md - Me = (48) 

for quarks. For both sets of particles /j,n c = —\i e . The 
baryon chemical potential can also be expressed in a sim- 
pler way: //jv b — M« for nucleons and Miv B = 2/J<d + M« 
for quarks. 
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conserved charge 


chemical potentials 




hadrons 


quarks 


N B 


lA B = (1 - Y p )n k n + Y P n* 


v% B = (2 - r p )/i5 + (i + y p K 




+Y P n k e + (Y L - Y p )fit 


+Y P n k + (Y L - Y p )tit 


Ya 


k fc fc 

I^Nr, ~ Mf — Me 


fc fc fc 

MiVrv ~ Mf — Me 


yk 
1 P 


M/V p — Mp — Mn — Ml/ + Me 


M/V p = Mil — Md — MtS + Me 


yk 


fc fc 

MiV r . — Mk 


MiV r , — Mf 



TABLE I: The local chemical potentials of the baryon number Nb, electric charge No, proton number N p and lepton number 
Nl in terms of the chemical potentials of the particles in one phases if the baryon number and all fractions are fixed locally. 
The second column is for a hadronic phase composed of neutrons, protons, electrons and neutrinos and the third column for a 
phase of up and down quarks and electrons and neutrinos. The results also apply for strange quark matter, with fid — For 
global baryon number conservation fi I Ng — hn b follows. If some of the fractions are conserved only globally and are no longer 
restricted by local constraints, the corresponding chemical potentials become equal, too: fj^ — (il/ . 



B. Role of Neutrinos 

Before we continue, we want to show why neutrinos 
do not have to be included in the construction of the 
non-neutrino part of an equation of state in terms of 
(T, tib, Yp), i-e. if the proton number is conserved. 

First we note, that for a single phase (T,riB,Y p ) are 
sufficient to fix all particle numbers but neutrinos. The 
non-neutrino part of the EoS would also not change, if 
neutrinos were not included as part of the thermody- 
namic system right from the beginning. At the same 
time, the neutrino contribution is also independent of 
the non- neutrino EoS: If the lepton fraction is conserved, 
i.e. if they are completely trapped, the neutrino density 
is directly specified by n v (T, fj, v ) = (Yl — Y p )ub- With- 
out lepton number conservation, /x„ = also directly sets 
the neutrino contribution. 

Neutrinos also do not have to be included in the con- 
struction of a mixed phase if we only allow global lep- 
ton number conservation. Then it always follows that 
= nl, 1 (the only quantum number of the neutrino is 
the lepton number so that eq. (pTTj) applies). In all cases, 
neutrinos can be treated as ideal Fermi-gases, so that 
their contribution to the EoS in the two phases is exactly 
the same. Thus it is also sufficient to study the pressure 
equilibrium without taking neutrinos into account. 

When the lepton fraction is conserved globally, the 
neutrino chemical potential can also always be taken out 
of the chemical equilibrium conditions. Because the neu- 
trino chemical potentials and all locally fixed fractions 
are equal in the two phases, terms proportional to Y\Mv 
cancel in the equilibrium conditions of other globally con- 
served charges. One finds that the same conditions for 
chemical equilibrium of the non-neutrino part are ob- 
tained, as if neutrinos were not present at all. We con- 
clude that neutrinos do not influence the phase equilib- 
rium between the two phases as they are distributed uni- 
formly over the entire system, if the lepton fraction is con- 
served globally. It is sufficient to calculate an equation of 
state for protons, neutrons and electrons (or quarks and 



electrons) in terms of (T, tib,Y p ) and this equation of 
state can be used for all possible conditions under which 
the neutrinos appear. 

Instead if Y p is not conserved, i.e. one has an EoS in 
terms of (T, tib) or (T, Ub,Yl) the neutrinos influence 
the rest of the matter via the condition for weak equilib- 
rium. The neutrino contribution needs to be taken into 
account for the evaluation of the non- neutrino EoS. Thus 
the non-neutrino EoS will also depend on whether lepton 
number is conserved or not, because the conditions for 
weak-equilibrium cqs. (|45p and (|46p . are different from 
those for beta-equilibrium, eqs. (|4"T|) and (|4"5|) . In the 
former case the lepton number is conserved and the neu- 
trinos arc determined by Yl . In the latter case neutrinos 
are completely untrapped and the neutrino contribution 
becomes trivial, [iv = 0. The same beta-equilibrium con- 
ditions eqs. (|4"T|) and (|4"5]) are found if neutrinos are not 
included in the thermodynamic description. In general 
the EoS depends on the conditions under which the neu- 
trinos appear, as soon as Y p is not conserved. 

Let us now discuss the meaning and some interest- 
ing consequences of the result /ij, = for matter with 
neutrinos but without lepton number conservation. If a 
particle i carries no conserved charges [a^ = V7) with 
cq. ((20)) one finds immediately that its chemical potential 
is zero. If neutrinos can be described as an ideal gas in 
equilibrium it follows that N v = 0, which means that the 
number of neutrinos equals the number of antineutrinos. 
Only if T = both contributions vanish. Non-conserved 
Yl would correspond to the situation when the neutrino 
mean free path is much larger than the size of the NS 
so that neutrinos can leave the NS freely. Thus the en- 
ergy of the system is not conserved, but can be carried 
away by neutrinos as long as they are abundant. As a 
logical consequence the NS would cool immediately to 
T = if weak reaction rates were fast enough (infinitely 
large emissivities) to allow to describe the neutrinos as an 
ideal gas as part of the thermodynamic description. In 
reality it takes some 10 5 years until the neutron star has 
cooled to a core temperature of ~ 10 keV and the photon 
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cooling era is reached. The neutrinos are far away from 
equilibrium, their emissivities have to be calculated, and 
the description of the cooling process requires detailed 
numerical simulations 12811. 



C. Mixed Phases 

If a mixed phase exists, it is crucial whether a charge 
is conserved globally or locally. In the following we will 
assume that each of the conserved charges is either con- 
served globally, or its fraction is conserved locally with 
equal values in the two phases, as before. Obviously, local 
constraints in the form = are equivalent to Y^ = 0, 
i.e. we can interpret local electric charge neutrality as 
a locally fixed charge fraction, too. Before we start to 
discuss all relevant combinations of locally and globally 
conserved charges, we will analyze the physical meaning 
of the different local constraints. 



Local charge neutrality 

Already in [3lj it was pointed out, that depending on 
the surface tension and the Debye-screening length, local 
charge neutrality might be the better approximation for 
bulk matter. If a large surface tension drives the system 
to sizes much larger than the Debye screening length, 
only a negligible small charged surface layer in the order 
of the Debye screening length remains and the bulk of the 
matter becomes locally charge neutral. Most calculations 
for the PT to quark matter indicate that this is indeed 
the case. 

If instead the surface tension is small so that the typical 
structures are smaller than the Debye screening length, 
global charge neutrality is the more reasonable assump- 
tion. This applies to the liquid-gas PT of nuclear matter, 
as the Debye screening length is large. Anyhow we in- 
clude the assumption of local charge neutrality for the 
nuclear PT in the following discussion, because it is in- 
structive and the corresponding equilibrium conditions 
can easily be devolved to other kind of phase transitions 
of nuclear matter. 

If one requires that both of the two phases have to be 
locally charge neutral, two different local electric charge 
chemical potentials appear in cq. ((20)) leading to differ- 
ent chemical potentials of all electrically charged parti- 
cles (electrons, quarks and protons) in the two phases. 
If one would do the full calculation including finite-size 
(surface and Coulomb) effects and without the constraint 
of local electric charge neutrality, the total chemical po- 
tential of charged particles would be shifted by the local 
electrical potential, e.g. p? v = [l 1 + eV 1 leading to full 
chemical equilibrium, p,^ = flp 1 , see [3l[. In the present 
article we are discussing infinite matter without Coulomb 
forces, so that the electrical potential cannot be deter- 
mined and the artificial inequality of the chemical poten- 
tials of charged particles cannot be resolved. 



Locally fixed Y v , Yl or ub 

There is no physical reason why the proton fraction, 
the lepton fraction or the baryon number density should 
be conserved locally or should be equal in the two phases, 
as there is no long range force between the two phases 
which is associated with these charges. This would imply 
that the readjustment of the local proton fraction, lep- 
ton fraction and/or baryon density does not take place. 
Therefore chemical equilibrium with respect to a lo- 
cally conserved charge is not established between the two 
phases. 

A non-vanishing locally fixed density (e.g. = rig = 
Ub) influences the condition for mechanical equilibrium 
so that pressure equilibrium is not obtained from the first 
and second law of thermodynamics any more. For these 
constraints a change of the subvolumes would imply a 
change of the local baryon numbers, too. Instead of pres- 
sure equilibrium only a combination of the local pressures 
and the local chemical potentials are equal in the two 
phases. Consequently, the pressure would change discon- 
tinuously at the phase transition. Thus, we will not use 
constraints of non- vanishing locally fixed densities. 

The likewise rather unphysical assumptions of locally 
fixed Y p or Yj, might be wanted because they allow to 
achieve a Maxwell construction of the mixed phase at 
the cost of only partial chemical equilibrium. Because of 
the additional conserved charges besides Nb local charge 
neutrality alone is not sufficient for that (as in the case of 
cold neutron stars) and at least one other of the conserved 
charges needs to be fixed locally. This is our motivation 
to investigate locally fixed Y p or Yt, and we will focus 
on this aspect when discussing different scenarios in the 
following subsections. 

Discussion of different cases 

In Tables Hfl fc IIIII all the relevant combinations of local 
and global conservation laws of the conserved charges for 
the construction of a mixed phase for isothermal mat- 
ter in supernovae, proto-NS and cold NS are listed. We 
assume that the densities and fractions are cither con- 
served globally, or locally in a form Yj = Y p n = Y p , 
Y[ = Y" = Yl, n z c = n# = n c = 0, n B = n'J = n B . 
The final equilibrium conditions are expressed in terms of 
the chemical potentials of the particles in the two phases, 
in Table HT1 for the liquid-gas phase transition and in Ta- 
ble [TTT] for the hadron-quark phase transition. 

1. Case I 

In case la, besides local charge neutrality the proton 
and lepton fractions are fixed locally and the system has 
only one globally conserved charge, the baryon number. 
The only internal degree of freedom of the two phases is 
the local baryon density. For this case table |T] expresses 
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case 


conserved de 
globally 


nsities / fractions 
locally 


equilibrium conditions 


construction of 
mixed phase 







n B, {Y p ), {Y L ),n c 


- 


direct 


la 


TLB 


Y p ,Y L ,nc 


(1 - YpK + Yp^'p + nl) + (Y L - Y p )ul = 
(1 - Y p )i£ + Ypbtf + fi") + {Yl - Y p )ul' 


Maxwell 


lb 


TLB 


Yl, nc 


u' n + Y Li ii = a'J + Ylu'J 


Maxwell 


Ic 


TLB 


Y p ,n c 


(1 - Y p )ni + Yp(p p + ui) = (1 - Y p )fi + Y P (rf + ni 1 ) 


Maxwell 


Id 


riB 


nc 


I _ II 

f*"n fl>n 


Maxwell 


Ila 


nB,Y L 


Y p ,n c 


(1 - Y p )a' n + Y 9 (j4 + u{) = (1 - Y p )i£ + Y 9 (j4* + ill 1 ) 

nt = ni 1 


Maxwell/Gibbs 


lib 


n B ,Y L 


nc 


Mn Mn 

A = m" 


Gibbs 


Ilia 


riB,Y p 


Y L , nc 


Mn + Y L nl = ni 1 + Ylu'J 

Mp Mn Mf > Me flp — Mn M^ ' Me 


Gibbs 


Illb 


riB,Y p 


nc 


Mri Mn 

Mp + Me = Mp 1 + A*" 


Gibbs 


IV 


n B ,Y L , Y p 


nc 


Mn Mn 

m£ = m" 

Mp + Me = m" + m" 


Gibbs 


V 


n B ,Y L ,Yp,nc 




Mri Mra 

M^ = m" 
Mp Mp 

Me = Me 1 


Gibbs 



TABLE II: Equilibrium conditions for the liquid-gas phase transition of nuclear matter for fixed temperature T, baryon number 
density ns and charge density nc- The lepton fraction Yl and proton fraction Y p are conserved in addition in some cases. 
These charge densities/fractions are fixed locally (with equal values in the two phases) or globally. If Y p is not conserved weak 
equilibrium (eqs. (|45[) & (|46[0 is established in both phases. If Yl is not conserved ul = mJ — is obtained, leading to the 
same equilibrium conditions as if neutrinos were not included in the thermodynamic system. 



the local chemical potentials belonging to the conserved 
charges Ng, Nc, N p and N L in terms of the chemical 
potentials of the particles in the phase. 

There is only one global chemical potential with the 
corresponding equilibrium condition: 

MWb = Mw B = M" B ■ (49) 

The condition which is shown in table |T] expresses that 
only combinations of particles can be exchanged which 
do not change the local proton and lepton fractions and 
are electrically charge neutral to maintain local charge 
neutrality. Thus in the liquid gas phase transition only a 
combination of 1 — Y p neutrons, Y p electrons and protons 
and Yl — Y p neutrinos can be exchanged freely between 
the two phases. In the hadron-quark phase transition 
only (1 + 5^,) up and (2 — Y p ) down quarks, Y p electrons 
and Yl —Y p neutrinos can be exchanged. 

As there is only one globally conserved charge Nb, 
the pressure is constant across the phase transition. The 
Maxwell construction can be used and all the results for 
G = 1 of section Hill can be applied. The chemical poten- 
tials of all particles are different in the two phases, as all 
particles contribute to the locally conserved fractions. 



We note that this case was already discussed in 
Ref. [47) ] . Equivalent equilibrium conditions were found 
and the same conclusions about the disappearance of the 
mixed phase in a compact star were drawn. 

In the following we will use the results of case la in 
table [J to derive the equilibrium conditions for all other 
cases. When a fraction is not conserved locally but only 
globally, the two local chemical potentials specify the new 
equilibrium conditions. In eq. (|26p it was deduced that 
by conserving Yy instead of Yj?, the two local chemi- 
cal potentials /i* become equal. If one of the charges 
is actually not conserved any more, this has to be seen 
as a global criterion. To minimize the free energy with 
respect to this charge, in eq. (|2"8|) it was derived that a 
non-conserved fraction leads to two new local constraints 
for the chemical potentials, = /j, 1 / = 0, which replace 
the two locally fixed fractions used before. 

In case lb Y p is no longer conserved. By setting 
M/v = the weak equilibrium conditions (|45[) respec- 
tively (|46|) of section HV Al are obtained which now have to 
be fulfilled in both phases. In Ic the conservation of Yl is 
lifted, leading to = [i]} = 0. If both fractions are not 
conserved any more as in case Id, the beta-equilibrium 
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n B , (Y p ), (Y L ),n c 


- 


direct 


la 


riB 


Y p ,Y L ,n c 


(1 - Y p )u n + Yp(fi p + /if) + (Y L - Yp)u» = 
(2 - Y v )a d + (1 + Y p )fi u + Y vt g + {Yl - Y p )a? 


Maxwell 


lb 


TIB 


Y L ,nc 


fin + Y L fi" = 2fi d + fi u + Y L fi? 


Maxwell 


Ic 


riB 


Y p ,nc 


(1 - Yp)a n + Y p (fi p + /if) = (2 - Y p )a d + (1 + Yp)a u + Y p (jg 


Maxwell 


Id 


riB 


nc 


fin = 2/ld + flu 


Maxwell 


Ila 


n B ,Y L 


Y p ,n c 


(1 - Y p )a n + Y p (fl p + /if) = (2 - Y p )a d + (1 + Yp)a u + Y p fi? 


Maxwell/Gibbs 


lib 


n B ,Y L 


n c 


fi n = 2 fid + fi u 
= M? 


Gibbs 


Ilia 


n B ,Y p 


Y L ,n c 


+ Yl/1v = 2^i d + /i„ + Yt/lJ 
/ip — /in — fl" + fie = flu — fid — fl? + fl? 


Gibbs 


Illb 


riB,Y p 


nc 


fin = 2fl d + /lu 

Mp + A*? = 2 /iu + fid + fi? 


Gibbs 


IV 


n B ,Y L , Y p 


nc 


fl n = 2^i d + /i u 

A 1 ^ = A 4 ? 

flp + fie = 2 fl u + fl d + fl? 


Gibbs 


V 


n B ,Y L ,Yp,nc 




fin = 2fl d + /lu 

Ai? = f*? 

flp = 2fl u + /id 

At? = /i? 


Gibbs 



TABLE III: As Table ILT1 but now for the hadron-quark phase transition. fi d — fi s is always valid. 



conditions (|47]l respectively (|48[) are obtained. These re- 
sults are independent of the other local or global conser- 
vation laws. Thus we do not need to discuss the non- 
conservation of lepton and/or proton fraction in the fol- 
lowing cases again. 

In cases la to Id the different conserved charges allow 
to rewrite the equilibrium condition fJ. J N = /ijv B m the 
simplified forms presented in tables |TT] & IIIII Case Id 
describes a cold, deleptonized neutron star. Only global 
baryon number conservation and local charge neutrality 
are considered. The well-known result of the equality of 
the neutron chemical potentials is found for the Maxwell 
construction of the liquid-gas phase transition. In all 
other Maxwell constructions the equality of the baryon 
chemical potential n% takes a different form and in- 
volves additional particles besides the neutrons. Because 
of the inequality of \i!y c in the two phases the chemi- 
cal potentials of the electrically charged particles always 
remain different in the two phases in all cases la to Id. 



2. Case II 

In case II lepton number and baryon number are con- 
served globally. The second equilibrium condition /i i Nl — 
fijj from the global conservation of lepton number leads 



to the equality of the neutrino chemical potentials. Neu- 
trinos are the only particles which can be exchanged be- 
tween the two phases, if the baryon number, the electric 
charge and the proton fraction are kept constant in both 
phases. 

Case Ila assumes locally fixed Y p and local charge neu- 
trality in addition. The same equilibrium condition as in 
case Ic in which Yj, is not conserved is obtained for the 
non-neutrino part of the EoS. Case Ic gives the same de- 
scription of the non-neutrino EoS as case Ila. If one does 
not include the neutrinos in the thermodynamic descrip- 
tion at all, the same condition as in Ic are found. Once 
more this shows explicitly that the non-neutrino EoS is 
independent of the neutrino contribution. As discussed 
before, the neutrinos can be calculated separately as long 
as Y p is conserved and Yl is not fixed locally. From this 
point of view, the non-neutrino EoSs of cases with fixed 
Y p and globally conserved Yl are equivalent to fixed Y p 
and non-conserved Yl- 

For the non- neutrino EoS G — 1 , and the mixed phase 
without neutrinos can be calculated with the Maxwell 
construction. The mixed phase would disappear under 
the influence of gravity in a hydrostatic configuration. 
But the inclusion of neutrinos leads to an interesting ef- 
fect on the mixed phase: The neutrino contribution is 
simply given by n„(T, /i„) = (Yl — Y p )ub- Thus for in- 
creasing baryon density also the neutrino density has to 
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increase. Therefore the neutrino pressure is not constant 
across the phase transition, which is in agreement with 
our general result for G = 2. If the pressure is used as 
the continuously varying variable and is changed strictly 
monotonic (e.g. in a compact star), a mixed phase ap- 
pears only because of the presence of neutrinos. 

It is very interesting to see, that all cases with local 
charge neutrality in which Y p or are conserved globally 
will lead to an extended mixed phase in a compact star. 
After the star has cooled to T = and became completely 
deleptonized, Y p and Yl are no longer conserved, and case 
Id will be reached. Consequently the mixed phase will 
disappear during the evolution of the star. We note that 
such a scenario has not been considered in the literature 
so far. 

If we compare case Ha to case la, we see that the same 
fractions and charges are conserved. In both cases only 
the Maxwell construction is needed, but in case Ha the 
unphysical assumption of a locally fixed lepton fraction 
is not used. Thus case Ha should be preferred over case 
la. 

We conclude that case Ha (or equivalently Ic for the 
non-neutrino EoS) is the most convenient scenario which 
leads for fixed Y p to the desired Maxwell construction of 
the system without neutrinos. All other cases with con- 
served proton fraction Y p involve more than one globally 
conserved charge for the non-neutrino EoS and the ex- 
plicit evaluation of phase equilibrium is necessary. Be- 
cause of the additional global conservation of the pro- 
ton fraction, even with local charge neutrality a simple 
Maxwell construction is not possible for matter in super- 
novae or proto-neutron stars. 

In case lib Y p is no longer conserved, so that the sep- 
aration of the neutrino EoS is not possible. The Gibbs 
construction has to be done with the inclusion of neutri- 
nos. Case lib is physically meaningful, as local charge 
neutrality is the only local conservation law. 



3. Case III 

The proton fraction is conserved globally in case III. 
The general equilibrium condition [i N = /ij^ shown by 
case Ilia in tables |TT] & [Hi] expresses that only a proton 
and an electron can be moved from one phase into the 
other, if at the same time a neutron and a neutrino are 
converted backwards. All other combinations of parti- 
cles would change the local baryon number, the electric 
charge or the lepton fraction. 

In case Ilia the neutrino EoS cannot be separated from 
the rest of the EoS, as the lepton fraction is conserved lo- 
cally so that the Gibbs construction has to be performed. 
If one is only interested to achieve the Maxwellian case 
(without further reasoning why locally fixed Yl instead 
of locally fixed Y p is assumed) , the easier case Ha can be 
applied instead. 

Local charge neutrality is the only local conservation 
law in case Illb, which is the proper description of a phase 



transition with sufficiently large surface tension between 
the two phases. The equilibrium conditions for the local 
baryon chemical potentials simplify compared to Ic, in 
which the proton fraction was fixed locally. 



4. Case IV 

The non-neutrino constraints of case IV are equivalent 
to those of case Illb. Case IV gives the correct descrip- 
tion of locally charge neutral matter with neutrinos as 
part of the thermodynamic description. 



5. Case V 

In case V local electric charge neutrality is not re- 
quired any more, so that all charges are conserved glob- 
ally. Eq. ([31) applies now for all particles in the liquid-gas 
phase transition, and the chemical potentials of all par- 
ticles become equal. For the hadron-quark phase tran- 
sition the hadronic chemical potentials directly fix the 
quark chemical potentials and vice versa. 

Yl and Y p are conserved globally in cases IV and V. 
Thus the equilibrium conditions remain the same if one 
or both of the fractions are no longer conserved. But ev- 
ery non-conserved fraction gives rise to two new stronger 
local constraints. They contain the information about 
the chemical equilibrium between the two phases with 
respect to this fraction used before, so that one of the 
equilibrium conditions in tables |TT] & IIIII becomes mean- 
ingless. 



6. Case 

Finally, we want to discuss case in which the phase 
transition is somewhat constructed by hand. In case 0, 
all state variables are fixed locally. No mixed phase has 
to be calculated, as also the baryon density is fixed lo- 
cally: n B = n 1 ^ = riB- If conserved at all, the three con- 
served fractions are fixed locally, too. In case there are 
no globally conserved charges so that no chemical equi- 
librium condition between the two phases is obtained. 
Thus, the two equations of state of the two phases can 
be calculated completely separately and the phase tran- 
sition point is then set by one freely selectable condition. 

Similarly to the Maxwell construction for the state 
variable p, the subvolumes of the two phases remain ar- 
bitrary at the transition point and are not constrained 
by the equilibrium conditions. Accordingly the extensive 
variables cannot be determined, too. On the other hand, 
the local intensive variables remain independent of the 
volumes of the two phases. Thus the chemical potentials 
of the locally conserved charges, the local pressure and 
the local temperature remain well-defined even without 
knowing the two subvolumes. The two phases can be 
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treated as independent single homogeneous phases with 
unknown volume. 

If pressure equilibrium is taken as the criterion for the 
determination of the phase transition point (which in this 
case is not a consequence of the first and second law of 
thermodynamics any more) the two phases can only co- 
exist at one special density n c £ ex if the other state vari- 
ables are kept constant. No extended mixed phase ap- 
pears, and all thermodynamic quantities but the pressure 
and the chemical potentials of electrons for conserved Y p 
and neutrinos for conserved Y p and Yl change discon- 
tinuously across the phase transition. If electrons and 
neutrinos can be treated as ideal gases /4 = follows 
from the conservation of Y p and nl = i^J from the con- 
servation of Y p and Yl ■ If instead of the baryon number 
density the pressure is used as the continuous variable 
no mixed phase forms, either. At the transition point 
pressure and thermal equilibrium are established, but at 
least chemical equilibrium of the baryons is not. Most 
importantly, therefore the thermodynamic potential, the 
free energy F = —pV + X^-^Mij behaves discontinu- 
ously. Because of the local constraints, the sum ^\ /i,iVj 
will not be equal in the two phases. If the pressure is 
used instead of the volume as one of the state variables, 
the free enthalpy will also behave discontinuously when 
the transition point is crossed. Thus, the thermodynamic 
potential can also not be used to determine which of the 
two phases exists before and which one after the phase 
transition. The second law of thermodynamics is vio- 
lated and thus some other additional criterion has to be 
considered for that. 

Besides pressure equilibrium every other possible coex- 
istence condition can be applied. The conclusions remain 
the same, and in general all the conditions will lead to a 
discontinuous thermodynamic potential. 

Instead one can use the continuity of the correspond- 
ing thermodynamic potential (the free energy F for 
(T, V, C^), the free enthalpy G for (T,p, C 7 )) as the phase 
coexistence criterion. At the point where the two poten- 
tials are equal the phase transition occurs. Before and 
after the phase transition only the phase with the lower 
potential is present. However, with this choice the pres- 
sure (and all other thermodynamic quantities but the 
state variables) will behave discontinuously. 

We note that case with fixed Y v would imply for 
the liquid-gas phase transition, that actually no phase 
transition occurs, because the same particles appear with 
equal densities in the two phases. Case is only relevant 
if some internal degrees of freedom remain which can be 
different in the two phases. 



D. EoS in (S/N B ,p,Y L ,n c ) 

Apart from numerical simulations, the EoS is usually 
calculated for the state variables (S/Nb,p, Yl,tic)- Such 
an EoS can be used directly for the description of typ- 
ical representative configurations of proto-ncutron stars 



which are characterized by a constant lepton to baryon 
ratio in a first approximation. 

Because the pressure is chosen to be the independent 
continuously changing variable, there exist only two pos- 
sibilities: Either the phase transition occurs only at one 
single pressure p coex , or a mixed phase of the two phases 
forms over an extended range in pressure. In the first 
case no mixed phase needs to be calculated. If ther- 
mal and mechanical equilibrium are required, the direct 
phase transition point can be found where the temper- 
atures of the two phases become equal T 1 = T 11 , with 
locally fixed = S/Ng = S/N B , p 1 = p n = p, 

Y[ = y/ J = Yl, rJfj = Uq = nc and (arbitrary) glob- 
ally conserved baryon number Nb which can be shared 
by the two phases. Pressure equilibrium is automatically 
given, as the pressure is one of the state variables which 
is set to equal values in the two phases. 

This case is the adiabatic equivalence to case of the 
isothermal phase transitions, in which all state variables 
but the volume were fixed locally. Here temperature equi- 
librium is chosen as the constraint which determines the 
coexistence point. Even though thermal equilibrium is 
enforced, the thermodynamic potential, which is the en- 
thalpy H = TS + ^ iVi/ii, will change discontinuously 
at the transition point. In contrast in the Maxwell con- 
struction of an isothermal phase transition the coexis- 
tence pressure is determined by the proper equilibrium 
conditions of the chemical potentials. We showed that 
this is equivalent to a continuous thermodynamic poten- 
tial. 

There exists nothing similar for an adiabatic process 
to the isothermal case Ha, which allows an easy con- 
struction of the non-neutrino mixed phase but leads to a 
second order phase transition with the inclusion of neu- 
trinos. As explained before, the non-neutrino part of the 
EoS depends on the neutrino fraction as soon as Y p is not 
conserved any more. Thus it is not possible to construct 
the non-neutrino mixed phase by means of an equal pres- 
sure Maxwell construction and treat the neutrinos inde- 
pendently. In the adiabatic case, in all scenarios (except 
for the direct phase transitions) the Gibbs construction 
has to be used and the mixed phase has to be calculated 
explicitly with the contribution of the neutrinos. 

The only physical meaningful local constraint is local 
charge neutrality, corresponding to the case of a large 
surface energy between the two phases. Globally con- 
served (S/Nb,Yl) is equivalent to globally conserved 
(S, Nb, Nl) for constant (arbitrary) baryon number Nb- 
Thus the conditions for chemical equilibrium are the 
same as in case lib. As Y p is not conserved, weak equi- 
librium is established in the two phases. Only neutrinos 
and neutrons can be exchanged independently between 
the two phases in the liquid-gas phase transition. For 
the hadron-quark phase transition two down and one up 
quark can be converted into a neutron. 

The possibility of a locally charge neutral extended 
mixed phase does not exist for cold dclcptonized neutron 
stars. So far such a kind of a mixed phase was not studied 
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in the literature. We will analyze the properties of this 
new mixed phase and examine its disappearance during 
the cooling of the neutron star in a forthcoming separate 
study (IJ. 

If local charge neutrality is lifted, no local charges ex- 
ist any more, and chemical equilibrium is expressed by 
case V, with weak equilibrium in the two phases in ad- 
dition. For the liquid-gas phase transition the chemical 
potentials of all particles are equal in the two phases. 

Because of the stronger constraint of local charge neu- 
trality, one can expect that the locally charge neutral 
mixed phase will extend over a smaller range in pressure. 
In both cases the appearance of the mixed phase assures 
that all thermodynamic variables are well defined and the 
first derivatives of the thermodynamic potential change 
continuously with the pressure. 



E. Strange Matter 

In strange quark matter, in addition to the up and 
down quarks, N s strange quarks arc part of the thermo- 
dynamic system. In principle, strange quarks carry the 
additional quantum number of strangeness. There exist 
two possibilities to handle this additional quantum num- 
ber: First, one can use the total strangeness of the system 
indeed as an additional conserved charge. If strangeness 
is not taken to be identical to zero, it is necessary to 
calculate the EoS for all possible strangeness fractions 
Ys = Ns/Nb- The strangeness chemical potential /xjv s 
would appear in addition to the chemical potentials of 
the other conserved charges. This approach was e.g. used 
in Ref. [12] to describe strangeness separation in heavy 
ion collisions. In this early work for the first time the 
phase transition to strange quark matter was described 
by two separate conserved charges, baryon number and 
strangeness, leading to the phase transition of a multi- 
component body. Here we will not discuss the scenario 
of conserved strangeness any further but will leave it for 
future discussion. 

Second, there exists a simpler and more commonly 
used description of strange matter, by assuming equi- 
librium with respect to strangeness changing reactions, 
i.e. strangeness is not conserved, so that /ijy s = 0. In 
this case the quantum numbers of the strange quark are 
identical to the ones of the down quark, so that [id = /j, s 
because of eq. (|30|) . Then all results presented for two- 
flavor quark matter in this article can also be applied to 
strange quark matter, with the baryon number given by 
Nb = 1/3(N U + Nd + N s ) and the electric charge number 
by N c = 1/3(2N U -N d - N s ) - N e . 

The only subtlety arises when Y p is conserved. First 
of all it is necessary to reconsider the meaning of Y p 
for strange matter. One possibility would be to inter- 
pret Y p as the net electric charge carried by baryons, 
N p = 1/3(2N U -Na- N s ), so that N. p = N e still gives 
charge neutrality. In combination with baryon number 
conservation, the conservation of Y p leads then to a fixed 



number of up quarks, but only the sum of down and 
strange quarks is fixed, i.e. reactions which change down 
into strange quarks are still in equilibrium. Thus it would 
be necessary to assume that these reactions happen on a 
much shorter timescale than reactions which change the 
number of up quarks (semi-leptonic reactions). 

The argumentation followed in Ref. [38[ is a differ- 
ent one. The EoS of state is calculated for fixed Y p . 
But within the application in a core-collapse supernova, 
quark matter appears only at such large densities and 
temperatures, that neutrinos are completely trapped and 
weak equilibrium is established. Hence Y p is actually 
not conserved but only Yl remains approximately con- 
stant. Within the numerical simulation for a given Yj, the 
proper Y p is determined. Under these conditions, weak 
equilibrium is a reasonable assumption and it is not nec- 
essary to argument with the different timescales of the 
different reactions. 

The whole discussion of this subsection applies also for 
hyperonic matter, i.e. hadronic matter with strangeness, 
and the conclusions are analog. 



V. SUMMARY & CONCLUSIONS 

In section[H]wc derived the chemical equilibrium condi- 
tions for the chemical potentials of the particles. We used 
a formulation, in which the two local parts Cq of the two 
phases k = I, II of the conserved charge Co = Cq + Cq 1 
and the local charge fractions Yf = C^/Cq of the con- 
served charges C T = + C^ 1 were chosen to be the 
independent degrees of freedom. Finally in eqs. pp|) . the 
chemical potentials of the particles were expressed by the 
chemical potentials of these degrees of freedom. We se- 
lected this special set of state variables because it can be 
used for all additional local constraints which are con- 
sidered in this article: locally fixed fractions with equal 
values in the two phases, or local charge densities which 
are zero (e.g. local charge neutrality). From this general 
formulation we continued with a particular set of local 
and global conservation laws and showed that the local 
chemical potentials of globally conserved charges with- 
out local constraints become equal, see eqs. ([25]) & ([25]) . 
I.e. if no local constraint exists for a charge, chemical 
equilibrium with respect to this charge is established. In 
eq. (|32p it became apparent, that hqCq (where /iq is the 
chemical potential of the charge Co) is equivalent to the 
Gibbs potential. 

With these results we showed in section Mil that the 
qualitative properties of the phase transition only depend 
on the number of globally conserved charges G and that 
they are independent of locally fixed fractions. In case 
there is only one globally conserved charge, G = 1, the 
simple Maxwell construction can be used for the calcu- 
lation of the mixed phase. The mixed phase will vanish 
in a static configuration under the influence of gravity, 
leading to a discontinuous phase transition with a con- 
tinuous thermodynamic potential. If there is more than 
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one globally conserved charge, G > 2, the Gibbs con- 
struction applies and every point inside the mixed phase 
has to be calculated explicitly. An extended mixed phase 
always forms, leading to a second order phase transition. 
The assumption of additional local constraints allows to 
achieve the Maxwell construction, if all the charge frac- 
tions are fixed locally. However, this does not work for 
an adiabatic phase transition, as the temperature has to 
vary with the density to keep the entropy constant. 

We applied these general findings to phase transitions 
in typical astrophysical systems in section [TV] We con- 
sidered the liquid-gas phase transition and the transition 
from hadronic to quark matter at high densities. It was 
shown that the results for up and down quark matter 
are the same as for strange quark matter with fid = ji s 
if equilibrium with respect to strangeness changing re- 
actions is assumed. First we focused on single phases, 
and derived that the equilibrium condition with a non- 
conserved fraction is given by setting the corresponding 
chemical potential to zero. The well known weak equi- 
librium conditions (eqs. (jl5|) & (pr6"|) ) are found if the 
proton fraction is not conserved, and the conditions for 
beta equilibrium (eqs. (|47|) & (|48|) ) if both the lepton 
and the proton fraction are not conserved. If the lepton 
fraction is not conserved = \i l J = 0, if neutrinos are 
included in the thermodynamic system. Furthermore it 
was shown that the neutrino EoS is independent of the 
rest of the matter as long as Y p is conserved and Yl is 
not fixed locally. 

For mixed phases it is crucial whether a local or a 
global conservation law is applied. We argued that local 
charge neutrality might give a better description of mixed 
phases, if the typical size of the structures is larger than 
the Debye screening length, due to a large surface tension 
between the two phases. Instead there is no physical mo- 
tivation for a locally fixed proton or lepton fraction, as no 
long range force is associated with these charges. How- 
ever, we used these assumptions to reduce the number of 
globally conserved charges. 

We classified the equilibrium conditions between the 
two phases for all relevant cases (see tables [H] & IIII[) . All 
of them were derived by applying the procedure described 
in section ITT] If a charge is conserved only globally its 
chemical potentials become equal in the two phases. For 
this we used the results from table U which expresses the 
chemical potentials of the conserved charges in terms of 
the chemical potentials of the particles, if all of the frac- 
tions are fixed locally. If Y p or Y p and Yl are no longer 
conserved, this is a global constraint leading to weak- or 
beta-equilibrium in both of the two phases. The case 
Id, in which Y p and Yl are not conserved, gives the well 
known Maxwell construction for cold-deleptonized neu- 
tron stars. 

Next, we discussed all the different cases. For fixed Y p 
case Ha was most interesting, as it allowed to use the 
Maxwell construction for the non-neutrino EoS with lo- 
cally fixed Y p as the only additional local constraint in 
addition to local charge neutrality. With the inclusion 



of neutrinos and globally fixed Yl, an extended mixed 
phase will form, causing a second order phase transition. 
The non-neutrino part of case Ha was equivalent to case 
Ic. Compared to case la, case Ha should be used prefer- 
ably, as the Maxwell construction is achieved without the 
assumption of a locally fixed lepton fraction. Case Ilia 
is also based on a locally fixed lepton fraction, but still 
requires the Gibbs construction. If the Maxwell construc- 
tion is not necessary, case Illb and IV show the equilib- 
rium conditions for local charge neutrality. In case V no 
local constraints are applied, resulting in the equality of 
the chemical potentials of identical particles in the two 
phases. 

For non-conserved Y p nothing analogous to case Ha 
does exist, as the non-neutrino part becomes dependent 
on the neutrino contribution. Either one assumes a lo- 
cally fixed lepton fraction (case lb) or the Gibbs con- 
struction has to be used. The only physical meaningful 
local constraint in the latter case is local charge neutrality 
(case lib). The same is true for an adiabatic phase tran- 
sition in which no Maxwell construction can be achieved 
at all. 

The possibility of a direct phase transition always ex- 
ists, by fixing all the state variables locally so that no 
mixed phase has to be calculated at all. However this 
is contradictory to the second law of thermodynamics as 
the thermodynamic potential changes discontinuously. 

It would be interesting to study the quantitative prop- 
erties of the mixed phase and the phase transition by ap- 
plying the general conditions found in the present article 
to specific equations of state for hadronic and quark mat- 
ter and to analyze the consequences for the proto-ncutron 
stars' evolution during the cooling process and its stabil- 
ity. Another promising investigation would be the impli- 
cations of the different assumptions for the mixed phase 
in the dynamical environment of a supernova, as done 
e.g. in 1 381 ] . 

The unphysical assumptions of locally fixed fractions 
have to be seen as a tool to construct a mixed phase 
in a thermodynamic consistent way, without the need of 
complicated calculations of the mixed phase. So far, the 
necessary conditions for the Maxwell construction were 
not derived explicitly in the literature. 

The assumption of local charge neutrality is instead 
physical meaningful and represents a simple way to sim- 
ulate strong Coulomb forces, without the need to include 
finite-size and Coulomb effects. Yet, local charge neutral 
mixed phases have not been analyzed in a lepton-rich, 
hot environment for the quark-hadron phase transition. 
Such locally charge neutral structures can grow almost 
arbitrary in size. One can expect significant changes of 
dynamical properties like the neutrino emissivities and 
opacities or the thermal conductivity. Furthermore local 
charge neutrality has an interesting effect on the evolu- 
tion of the mixed phase itself. As long as Y p or Yl are 
conserved at least two globally conserved charges exist 
and a mixed phase is present. After the deleptoniza- 
tion, finally only Nb will remain as a globally conserved 
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charge, and when the star becomes isothermal the mixed 
phase will disappear. The authors are currently prepar- 
ing a separate study [5l| in which these effects will be 
analyzed for some specific EoSs. 

There is another application of the results in this arti- 
cle, namely the description of the nucleation of the quark 
phase. As the initially formed bubble will have a small 
size, surface effects become crucial. The deconfinement 
process is mediated only by strong interactions. Since the 
timescales for any weak reactions are much longer, one 
can assume that this leads to flavor conservation during 
the deconfinement transition. Depending on the size of 
the relevant fluctuations in the hadronic phase one might 
choose one of the different cases presented in table Mil to 
calculate the threshold density of the nucleation. E.g. in 
the equilibrium conditions denoted by case la in the 
present article were derived with the motivation to de- 
termine the phase transition point under the assumption 
of flavor conservation. The authors assume locally fixed 
Y p and Yl and local charge neutrality. However, in this 
calculation chemical equilibrium with respect to baryon 
number is still established, leading to different baryon 
baryon densities in the two phases, although the proton 
fraction remains equal in the two phases. One might con- 
sider the even stronger constraint, that initially also the 
baryon density has to remain constant inside the bubble 
which is nucleated. After the quarks have been decon- 
fined, weak reactions will take place and the system will 
develop to a state described by case IV for local charge 
neutrality or by case V for global charge neutrality. La- 
tent heat will be released during this non-equilibrium 
process. 

The present work may also be relevant for the descrip- 
tion of the QCD phase transition in relativistic heavy ion 
collisions, explored in experiments at e.g. CERN, RHIC 
or FAIR. If the relativistic expansion of the fireball in 



the central collision zone is characterized by a constant 
entropy and baryon number, from our considerations we 
can conclude that a mixed phase exists over an extended 
range in pressure and thus over a finite period in time, 
too. The phase transition will be continuous, if ther- 
modynamic equilibrium is reached at all stages of the 
expansion. The description of matter in the QCD phase 
transition as a multi-component system allows different 
choices of the equilibrium conditions, as presented here 
in the case of compact stars. E.g. in Ref. [53| in addi- 
tion to global conservation of baryon number the global 
conservation of isospin was considered. The recent inves- 
tigation [H[ indicates, that also strangeness fluctuations 
and separation could lead to interesting effects. However, 
in contrast to the case of compact stars, the typical hy- 
drocells in a heavy ion collision arc usually very small of 
order 0.1 fm. Thus it might be questionable to use an 
EoS which includes a mixed phase for the hydrodynamic 
description. We note that only such globally conserved 
charges are relevant degrees of freedom, which are actu- 
ally accessible by the system. 
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